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Abstract. Energy spectrum and thermodynamics of the pseudospin-electron model introduced at the 
consideration of the anharmonicity effects in high-T c superconductors are investigated in the dynamical 
mean field approximation (d — oo limit). In the limit of zero electron correlation U — > this model is 
analytically exactly soluble within this approach: in the fx = const regime the first order phase transition 
with the jump of the pseudospin mean value (S z ) and reconstruction of the electron spectrum can realize, 
while in the n = const regime the phase separation in electron subsystem can take place for certain values 
of the model parameters. On the basis of the obtained results the applicability of the approximate schemes 
previously used for the investigation of the pseudospin-electron model are discussed. 

PACS. 71.10.Fd Lattice fermion models (Hubbard model, etc.) - 71.38.+i Polarons and electron-phonon 
interactions - 77.80.Bh Phase transitions and Curie point - 63.20.Ry Anharmonic lattice modes 



1 Introduction 

The theoretical investigation of the strongly correlated 
electron systems is an enduring subject of interest in con- 
densed matter physics especially during last ten years after 
the discovery of high-T c superconductivity. Recent stud- 
ies of strongly correlated Hubbard type models elucidate 
some important features of high-T c , namely the G?-wave 
pairing and the role of antiferromagnetic fluctuations [EJ . 
However some features of the cuprates are not well un- 
derstood (e.g. the existence of an "optimal" doping, the 
effect of charge fluctuations, strong electron-phonon in- 
teraction, lattice dynamics and instabilities of ferroelec- 
' trie type). Thus the extensions of the Hubbard model by 
, the van Hove scenario concept, the incorporation of the 
local electron-phonon interaction (the Hubbard-Holstein 
model for harmonic phonons and the pseudospin-electron 
model (2) for anharmonic ones), generalization to a two 
or three band model, the inclusion of intersite electron 
interaction, etc. are under consideration. 

Within these models the pseudospin-electron one @ 
in a simplest way includes the interaction of correlated 
electrons with some local lattice excitations described by 
pseudospins (e.g. anharmonic vibrations of apex oxygen 
in YBaCuO type HTSC's), and shows the possibility of 
dipole (pseudospin) and charge density instabilities ||,^| 
and phase separation || due to the effective retarded inter- 
action between pseudospins via conducting electrons. All 
these results were obtained within the generalized random 
field approximation (GRPA) Q which is a realization of 
the appropriate perturbation theory for correlation func- 
tions in the case of strong coupling (U 3> t) and corre- 
sponds to the mean field approximation in calculation of 
mean values. There are no good criteria of its applicability 



and it is supposed that GRPA gives correct description in 
the case of large dimensionality of local (site) states. 

In recent years the essential achievements of the the- 
ory of strong correlated electron systems are connected 
with the development of the dynamical mean field theory 
(DMFT) proposed by Metzner and Vollhardt @ for Hub- 
bard model (see also || and references therein). DMFT 
is a nonperturbative scheme which allows to project 
Hubbard model on the single impurity Anderson model 
and is exact in the limit of infinite space dimensions 
(d = 00). Moreover, some class of models (e.g. Falicov- 
Kimball model §) can be studied almost analytically 
within DMFT. 

Here we apply DMFT to the investigation of pseudo- 
spin-electron model in the limit of zero electron correla- 
tion (U — 0) which can be treated analytically. 



2 Perturbation theory in terms of electron 
transfer 

The Hamiltonian of pseudospin-electron model in the ab- 
sence of electron correlations can be written in the form: 



1 j" 



where 



Hi = gSf y^ur - n) j n lcr - hS\ 



(1) 



(2) 



is single-site Hamiltonian, and includes local interaction 
of conducting electrons with pseudospins placed in longi- 
tudinal field h (asymmetry parameter of anharmonic po- 
tential) . 
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In general, one-electron Green's function G a (<jJ n , k) 
G%(j r') = (Ta itT (T)a] a (r')a((3)) o / (a((3)) (3) 

a(P) = T exp | - J\t g t y at(r)a JCT (r)| (4) 
satisfies Larkin's equation 

G%{T-r') = S^T-T / ) + S! l (T-T")t lm G^(r"-T'), (5) 

where summation (integration) over repeated indices is 
supposed. The formal solution of eq. (^) can be written in 
the form 

G a (u n ,k) = — — (6) 

(o; n ,k)-i k 

and the task is to calculate the irreducible according to 
Larkin parts S a (u) n ,'k). 

It is convenient to introduce projective operators on 
pseudospin states 



~±Sf, (P±) 2 = P?, P t + P t -=0 (7) 



and by substitution = a, P^ = 1 — Cj Hamiltonian 
([!]) , (0) can be transformed into the Hamiltonian of binary 
alloy. On the other hand, if we keep in ([!]), (||) only elec- 
trons with one orientation of spin by removing the sum 
over spin indices and putting a =| and consider electrons 
with a =\ as localised Pf = n$j_, P~ = 1 — we get the 
Hamiltonian of the Falicov-Kimball model where h plays a 
role of the chemical potential for the localized .[-electrons. 
As a rule the common chemical potential is introduced 
for both electron subsystems but the case of two chemical 
potentials was also considered (see, e.g. (l(J) and the first 
consideration of Falicov-Kimball model within DMFT was 
done by Brandt and Mielsch |J. 

Formally diagrammatic series for the irreducible part 
3?j (u> n ) are the same for all these models 



SS(<0») 




(8) 



and includes both single-site and intersite contributions. 
Here, arrows indicate electron propagators = 

p ± 

iuj +'/jt- ' n ^ e su bspaces projected on the pseudospin 
states of site i and ovals represent semi-invariant averag- 
ing of projection operators. 

The main difference between these models is in the 
way how an averaging procedure over projection opera- 
tors is performed (thermal statistical averaging in the case 
of pseudospin-electron and Falicov-Kimball models and 
configurational averaging for binary alloy) and how self- 
consistency is achieved (fixed value of longitudinal field h 



for pseudospin-electron model, fixed value of the compo- 
nent concentration c for binary alloy and fixed value of 
the electron concentration — total or for both electron 
subsystems — for Falicov-Kimball model). 



3 The limit of large dimensions (d — > oo) 

oo) one should scale 



In the case of high dimensions (d 
hopping integral 

i » hi. 

ij Vd 



(9) 



in order to obtain finite density of states (the Gaussian 
one for d = oo hypercubic lattice p(e) = ^^ e -6 > w 
and semi-elliptic d.o.s. for d = oo Bethe lattice p(e) = 
-^z^fW 2 — e 2 (|]). Due to such scaling only single-site 
contributions survive in the expression for irreducible 
parts S a 

- t') = b l3 E a {T - t'), E a (cj n , k) = 3 a {u n ) (10) 

and such site-diagonal function, as it was shown by 
Brandt and Mielsch || , can be calculated by mapping the 
infinite-dimensional lattice problem on the atomic model 



-/3H 



-0H ei! 



-0H O 



(11) 



x T exp | — 



dr I dr 



'E 



J a {r - T')4(T)a a (r') 



with auxiliary Kadanoff-Baym field J a (r — r') |Tl| which 
has to be self consistently determined from the condition 
that the same function S a defines Green's functions for 
lattice (p) and atomic limit 



5 a (u> n ) - Jo-(Wn) 



(12) 



"Dynamical" mean field J CT (r — r') describes the hop- 
ping (transfer) of electron from atom into environment at 
moment r, propagation in environment without stray into 
atom until moment r'. Connection between these "dynam- 
ical" mean field of atomic problem and Green's function of 
the lattice can be obtained using standard CPA approach 



where 



dt- 



P(t) 



(13) 



(14) 



is a single-site Green's function both for atomic limit and 
lattice. Here summation over wave vector was changed by 
the integration with the density of states p(t). 

In order to complete our self-consistent set of equa- 
tions we should find expression for Green's function in 
the atomic limit (|l2|). Due to the properties of the projec- 
tion operators (0) one can rewrite Hamiltonian of atomic 
problem (111]) in the form 



-PH<a — p+ e -0H l+) 



(15) 
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and our space of states splits into two independent sub- 
spaces hence all projection operators (0) act at the same 
site and in any order of the perturbation theory expansion 
all projection operators can be replaced by their product 
result and there are no necessity to make semi-invariant 
expansions. 

Single-electron Green's function is a sum of Green's 
functions in subspaces and is equal 



+ 



iu) n + fi— Ja(w n ) - f 

iuj n J a {uJ n ) + § 



(16) 



Partition functions in subspaces are 
Z ± = Spe-^±=e ± ^-^ (17) 



n(i+e-^«)ni- 



iuj n + \x =F f 

and presents the partition functions of the non-interacting 
fcrmions with frequency dependent hopping placed in the 
external field formed by pseudospin. 

Pseudospin mean value is determined by equation 



(s z ) = l Z+ Z 



2Z+ + Z- 
\^±{flh-(Q + [{S')]-Q-[{S')])) 



(18) 



which is an analogue of the well known equation of state 
for Ising model in mean-field approximation: (S z ) = 
g tanh ^ (ft, + Jq (S z )). It should be noted that in the case 
of Lorentzian density of states p(e) = ^y$z+gi\ , which is 
frequently used in some applications of DMFT, one can 
easily obtain a simple result J a (uj n ) = i W M , quantities 
Q ± do not depend on (S z ) and equation ( |l9| ) transforms 
into an expression for (S z ) that indicates the sensitivity 
of the equation of state to the shape of d.o.s. 

Electron concentration mean value is determined by 



(19) 



and the functional of thermodynamic potential can be de- 
rived in the same way as it was done in § for Falicov- 
Kimball model 

^ = %) - \ E { lnG S 0) K) m G a (uj n , k) 

^ Tl.fT ^ lr 



where 



tt(«) = --gln(Z+ + Z_) 



(20) 
(21) 



is a thermodynamic potential for atomic problem. 

Below, all calculations will be performed for semi- 
elliptic density of states when the auxiliary field is de- 
termined by the simple cubic equation 



W 2 



(P+) 



iw n + fj, - J a (u n ) 



(22) 
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Fig. 1. Electron bands boundaries (semi-elliptic d.o.s., W = 
0.4, (S z ) = 0.2). 



+ fi - J ff (w„) + 



In a usual way we perform analytical continuation on 
real axis (iu) n — > w — iS) and only solutions of (p2|) with 
5mJ (7 (ci;) > must be considered. Band boundaries are 
determined from the condition SmJ 5 (u) — > and in Fig. 1 
their dependence on coupling constant g are presented. 
One can see that there exists critical value of coupling 
constant g ~ W when a gap in spectrum appears. It 
should be noted that within GRPA as well as in other 
approaches where single-electron Green's function is cal- 
culated in Hubbard-I approximation, when we keep only 
the first term of the single-site contribution in the expres- 
sion for the irreducible part (|^), these gap in spectrum 
always exists. 

In the case of strong coupling (g 3> W) an analytical 
solutions can be obtained 



(23) 



for upper and lower subbands, respectively, and one can 

see that subbands halfwidth is equal to W | ± (S z ) 

whereas in Hubbard-I approximation it is W (4 i (S z ))- 
This result clearly shows that even for the case of strong 
coupling when subbands are well separated and one of 
them become narrow ((S z ) — > ±4) Hubbard-I approxima- 
tion is unsufneient and can not be derived from the exact 
solution in any way, e.g. due to the subbands halfwidth 
square root dependence on the localized states occupance 
«P*>-0). 

Presented above expressions were obtained for the 
fixed value of the chemical potential /i when stable states 
are determined from the minimum of the thermodynami- 
cal potential (|20|). This regime fi = const corresponds to 
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the case when the charge redistribution between conduct- 
ing sheets CuC>2 and other structural elements (charge 
reservoir, e.g. nonstoichiomctric in oxygen CuO chains in 
YBaCuO type structures) which fix the value of the chem- 
ical potential is allowed. On the other hand, in the regime 
of the fixed electron concentration value one should solve 
equation for chemical potential n = (n) (|l9|) and stable 
states are determined by the minimum of the free energy 
F = Q + im. 



4 Results and discussion 

Integrals in Eqs. ( |l7| ) and (^0|) can be calculated analyti- 
cally for states with (S z ) = ±| at zero temperature and 
corresponding phase diagrams \x—h which indicate stabil- 
ity regions for these states are shown in Fig. |^a and b for 
g > W and g < W, respectively. One can see two regions 
of n and h values where the states with (S z ) = ±i coex- 
ists. In the vicinity of these regions the phase transitions 
of first order with the change of the longitudinal field h 
and/or chemical potential fj, take place (see Fig. ||) and 
they are shown by thick lines on phase diagrams (Fig. |J) . 

There are no any specific behaviour when chemi- 
cal potential is placed out of bands. If chemical poten- 
tial is placed in upper subband the graphs presented in 
Fig. [| transform according to the internal symmetry of 
the Hamiltonian: 



-fj,, 



2g — h, n 



n, 



S z 



-S z . (24) 



With the temperature increase the region of the phase 
coexistence narrows and the corresponding phase diagram 
T c — h is shown in Fig. ||. One can see that with respect to 
Ising model the phase coexistence curve is shifted in field 
and distorted from the vertical line and hence the possibil- 
ity of the first order phase transition with the temperature 
change exists in pseudospin-electron model for the narrow 
range of h values. 

As it was mentioned above, the band structure is de- 
termined by the pseudospin mean value and its change is 
accompanied by the corresponding changes of the electron 
concentration and for the (//, h) values fixed on the first 
order phase transition line there are three solutions for 
electron concentration one of which is unstable. 

In the case of the fixed value of the electron concentra- 
tion value (regime n — const) this first order phase tran- 
sition transforms into the phase separation. One can see 
regions with dfi/dn < 0, which correspond to this effect 
in electron subsystem, on the concentration dependencies 
(Figs. |and |a). 

The corresponding dependencies of free energy F = 
Q + an are given in Fig. |b. In the phase separated region 
free energy deflects up and concentration values at binodal 
points are determined by the tangent line touch points or 
from the chemical potential dependencies (Fig. ga) using 
Maxwell construction. Resulting phase diagram T — n is 
shown in Fig. [?]. 
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Fig. 2. Phase diagram n — h. Dashed and thin solid lines 
surround regions with S z = respectively. Thick solid line 
indicate the first order phase transition points, a) g = 1, W = 
0.2; b)g = l,W = 0.7. 



For the first time the possibility of phase separation in 
pseudospin-electron model was marked in [^| where it was 
obtained within GRPA in the limit of strong correlation 
U — > oo. Here it is observed for the opposite case of U = 0. 

The problem of phase separation in strongly corre- 
lated systems is not new (see and references therein). 
It was shown for Hubbard and t — J models that 
for some parameter values system separates into hole-rich 
and hole-poor regions with paramagnetic and antiferro- 
magnetic orders, respectively, and long-range interaction 
between these charged regions is considered as an origin of 
the appearance of stripe structures. In Rcf. |l(| the phase 
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Fig. 3. Field dependencies of (S z ) (a) and thermodynamical 
potential (b) for [i — const regime when chemical potential 
is placed in the lower subband fi — —0.37 (W = 0.2, g = 1, 
T = 0). 



segregation for some parameter values was reported for the 
annealed binary alloy with diagonal disorder described by 
Falicov-Kimball model. In our case of pseudospin-electron 
model without electron correlations system separates into 
regions with different values of electron concentration and 
pseudospin mean value and electron spectrum contains 
both wide empty electron band and occupied localized 
states of the regions with n ~ as well as partially filled 
wide electron band and empty localized states of the re- 
gions with n ~ 1 (see Fig. ||) the weights of which are 
determined by the electron concentration. Such type local- 
ized states (polarons) results from the strong electron-out 
of plane apical oxygen vibrations coupling (g > W) in the 
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Fig. 4. Phase diagram T c — h: solid and dashed lines indicate 
the first order phase transition line and boundaries of the phase 
stability, respectively (g = 1, W = 0.2, fi = —0.5) 
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Fig. 5. Dependence of the chemical potential p and electron 
bands boundaries (dashed lines) on the electron concentration 
n(T = 0.001, g = 1, W = 0.2, h = 0.1). 



case of YBaCuO-type structures and it is supposed that 
the hopping between them gives significant contribution 
in the carrier relaxation observed by the resonant Raman 
spectroscopy Q. 

It should be noted that in the case of spinless fermions 
Hamiltonian (|l|) can be applied for the description of the 
oxygen vacancies subsystem in high-T c superconductors, 
which can be treated as quasiequilibrium, and it is known 
that their interaction with some relaxation type lattice 
mode leads to the phase separation and appearance of 
superstructures and stripes UM . 
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Fig. 6. Dependence of the chemical potential \i (a) and de- 
viation of free energy from linear dependence AF — F(n) — 
j F(2) — (l — § ) .F(O) (b) on the electron concentration n for 
different temperatures T (g = 1, W = 0.2, h = 0.1). 



In this paper we investigated the possible phase transi- 
tions in pseudospin-electron model within DMFT without 
creation of super structures (k = 0) and the phase dia- 
grams presented in Figs. |^ and ^ concern only this case. In 
order to detect instabilities associated with a specific wave 
vectors one should calculate response functions which will 
be the subject of the further investigations. 
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Fig. 7. Phase diagram T — n for phase separated state: solid 
line — binodal, dashed line — spinodal (g = 1, W = 0.2, 
h = 0.1). 
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